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Abstract 

In a previous study (Matsuno Y J. Phys. A: Math. Theor. 45 (2012) 23202), we have 
developed a systematic method for obtaining the bright soliton solutions of the Fokas- 
Lenells derivative nonlinear Schrodinger equation (FL equation shortly) under vanishing 
boundary condition. In this paper, we apply the method to the FL equation with non- 
vanishing boundary condition. In particular, we deal with a more sophisticated problem 
on the dark soliton solutions with a plane wave boundary condition. We first derive the 
novel system of bilinear equations which is reduced from the FL equation through a de- 
pendent variable transformation and then construct the general dark iV-soliton solution of 
the system, where N is an arbitrary positive integer. In the process, a trilinear equation 
derived from the system of bilinear equations plays an important role. As a byproduct, 
this equation gives the dark iV-soliton solution of the derivative nonlinear Schrodinger 
equation on the background of a plane wave. We then investigate the properties of the 
one-soliton solutions in detail, showing that both the dark and bright solitons appear on 
the nonzero background which reduce to algebraic solitons in specific limits. Last, we 
perform the asymptotic analysis of the two- and iV-soliton solutions for large time and 
clarify their structure and dynamics. 
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1. Introduction 



The Fokas-Lenells derivative nonlinear Schrodinger (NLS) equation (FL equation shortly) 
is a completely integrable nonlinear partial differential equation (PDE) which has been 
derived as an integrable generalization of the NLS equation using bi-Hamiltonian methods 
[1]. In the context of nonlinear optics, the FL equation models the propagation of non- 
linear light pulses in monomode optical fibers when certain higher-order nonlinear effects 
are taken into account [2]. We employ the following equation which can be derived from 
its original version by a simple change of variables combined with a gauge transformation 
[2]: 

u x t = u — 2i\u\ 2 u x . (1.1) 

Here, u = u(x, t) is a complex-valued function of x and t, and subscripts x and t appended 
to u denote partial differentiations. The complete integrability of the FL equation has been 
demonstrated by means of the inverse scattering transform (1ST) method [3]. Especially, 
a Lax pair and a few conservation laws associated with it have been obtained explicitly 
using the bi-Hamiltonian structure and the multisoliton solutions have been derived by 
applying the dressing method [4]. Another remarkable feature of the FL equation is that 
it is the first negative flow of the integrable hierarchy of the derivative NLS equation [2, 
5]. 

In a previous study [6] which is referred to as I hereafter, the two different expressions 
of the bright iV-soliton solution of the FL equation have been obtained by a direct method 
which does not recourse to the 1ST and their properties have been explored in detail. Here, 
we construct the dark 7V-soliton solution of the FL equation on the background of a plane 
wave. Explicitly, we consider the boundary condition 

u — > pexp {i (kx — cot + ( - ± - ) ) } , x — > ±oo, (1.2) 

where p{> 0) and k are real constants representing the amplitude and wavenumber, 
respectively, <fi^ are real phase constants and the angular frequency u = oo{k) obeys the 
dispersion relation oj = 1/k + 2p 2 . Note that the plane wave given in (1.2) is an exact 
solution of the FL equation. As will be discussed later, the possible values of k must be 
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restricted to assure the existence of the soliton solutions. A similar problem to that posed 
in this paper has been studied recently and an explicit formula for the dark iV-soliton 
solution have been presented by an ingenious application of the Backlund transformation 
between solutions of the FL equation and the Ablowitz-Ladik hierarchy [7]. Nevertheless, 
the detailed analysis of the soliton solutions has not been undertaken as yet. 

An exact method of solution employed here which is sometimes called the direct 
method [8] or the bilinear transformation method [9] is a powerful tool for analyzing 
soliton equations and differs from the method used in [7]. Once the equation under con- 
sideration is transformed to a system of bilinear equations, the standard technique in the 
bilinear formalism is applied to obtain soliton solutions. A novel feature of the bilineariza- 
tion of the FL equation is that one of the bilinear equations can be replaced by a trilinear 
equation, as already demonstrated in I. The same situation happens in the current dark 
soliton problem. However, the resulting trilinear equation will be used essentially in the 
process of performing the proof of the dark iV-soliton solution. 

This paper is organized as follows. In section 2, we bilinearize the FL equation under 
the boundary condition (1.2). We then show that one of the resulting bilinear equations 
can be replaced by a trilinear equation. In section 3, we present the dark iV-soliton so- 
lution of the bilinear equations. It has a simple structure expressed in terms of certain 
determinants. Subsequently, we perform the proof of the dark iV-soliton solution using an 
elementary theory of determinants in which Jacobi's identity will play a central role. As 
already noted, the proof of the trilinear equation turns out to be a core in the analysis. 
In accordance with the relation between the FL equation and the derivative NLS equa- 
tion at the level of the Lax representation, we also demonstrate that the dark iV-soliton 
solution obtained here yields the dark iV-soliton solution of the derivative NLS equation 
by replacing simply the time dependence of the solution. As in the case of the defocusing 
NLS equation subjected to nonvanishing boundary conditions, it is necessary for the ex- 
istence of dark solitons that the asymptotic state given by (1.2) must be stable. Hence, 
we perform the linear stability analysis of the plane wave solution (1.2) and provide a cri- 
terion for the stability. In section 4, we first investigate the properties of the one-soliton 
solution in detail. We find that depending on the sign of k and that of the real part of the 
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complex amplitude parameter, the solution can be classified into two types, i.e., the dark 
and bright solitons. The latter soliton may be termed " anti-dark soliton" since the back- 
ground field is nonzero. However, we use a term "bright soliton" throughout the paper. 
We demonstrate that regardless the sign of k, the bright soliton has a limiting profile of 
algebraic type (or an algebraic bright soliton) whereas an algebraic dark soliton appears 
only if k < 0. We then analyze the asymptotic behavior of the two-soliton solution and 
derive the explicit formulas for the phase shift in terms of the amplitude parameters of 
solitons. In particular, we address the interaction between a dark soliton and a bright 
soliton as well as that of two dark solitons. Last, the similar asymptotic analysis to that 
of the two-soliton solution is performed for the general dark iV-soliton solution. Section 
5 is devoted to concluding remarks. 

2. Exact method of solution 

In this section, we develop a direct method of solution for constructing dark soliton 
solutions of the FL equation (1.1) under the boundary condition (1.2). In particular, we 
show that it can be transformed to a system of bilinear equations by introducing the same 
type of the dependent variable transformation as that employed in I for the bilinearization 
of the FL equation under vanishing boundary condition. We also demonstrate that this 
system yields a trilinear equation which will play a crucial role in our analysis. 

2.1. Bilinearization 

The bilinearization of the FL equation (1.1) is established by the following proposition: 
Proposition 2.1. By means of the dependent variable transformation 



with ou — 1/k + 2p 2 , equation (1.1) can be decoupled into the following system of bilinear 
equations for the tau functions f and g 



u = pe 



(2.1) 



D t f ■ f* - v 2 (gg* - ffl = o, 



(2.2) 



D x D t f ■ f* - ip 2 D x g ■ g* + ip 2 D x f ■ f* + 2np 2 (gg* - //*) = 0, 



(2.3) 
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D x D t g ■ f + mD t g ■ f - iuD x g ■ f = 0. (2.4) 

Here, f = f(x,t) and g = g(x,t) are complex-valued functions of x andt, and the asterisk 
appended to f and g denotes complex conjugate and the bilinear operators D x and D t are 
defined by 

where m and n are nonnegative integers. 

Proof. Substituting (2.1) into (1.1) and rewriting the resultant equation in terms of the 
bilinear operators, equation (1.1) can be rewritten as 

^{D x D t g ■ f + mD t g ■ f - iuD x g ■ f) 
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{f*D x DJ ■ f - 2np 2 fr ~ 2ip 2 g*(gJ - gf x + mfg)} = 0. (2.6) 



Pf 

Inserting the identity 

f*D x D t f ■ f = fD x D t f ■ f* ~ 2f x DJ ■ f* + /(A/ • Dx, (2-7) 

which can be verified by direct calculation, into the second term on the left-hand side of 
(2.6), one modifies it in the form 

j^(D x D t g ■ f + \KD t g ■ f - iuD x g ■ f) 
-Jp [f{D x D t f ■ f* - ip 2 D x g ■ g* + \p 2 D x f • /* + ^p 2 (99* ~ //*)} 

-2/4 a/ • r - ip 2 (gg* - //*)} + f{DJ ■ f* - ip 2 (gg* - ff*)} x ] = o. (2.8) 

By virtue of equations (2.2)-(2.4), the left-hand side of (2.8) vanishes identically. □ 
It follows from (2.1) and (2.2) that 

H 2 = p 2 + i|ln^. (2.9) 
The above formula gives the modulus of u in terms of the tau function /. 



5 



2.2. Trilinear equation 

Proposition 2.2. The trilinear equation for f and g 

f* \g*tf - {fx - i«/)<& - i + p 2 ) {gj - g f x ) } = f t *{g x f - g fx + infg), (2.10) 

is a consequence of the bilinear equations (2.2)-(2.4). 

Proof. By direct calculation, one can show the following trilinear identity among the tau 
functions / and g: 

f* ^9xtf ~ {fx - iKf)g t - i Q + p 2 ^j {g x f - gf x )^ - f* t {g x f - gf x + infg) 



= f*{D x D t g ■ f + iKD t g ■ f - \uD x g ■ f) 

~{DJ ■ f* - ip 2 {99* - //*) } x + {DxD t f ■ f* - ip 2 D x g ■ g* + \p 2 D x f ■ f* - 2«£>*/ • /* 

+ g x {D t f.f*-ip 2 {gg*-ff*)}. (2.11) 

Replacing a term 2inD t f ■ f* on the right-hand side of (2.11) by (2.2), the right-hand side 
becomes zero by (2.2)-(2.4). This yields (2.10). □ 

In view of proposition 2.2, the proof of the dark iV-soliton solution is completed if 
one can prove any three equations among the three bilinear equations (2.2)-(2.4) and a 
trilinear (2.10). We will see later in section 3 that the proof of (2.4) is not easy to perform 
and hence we prove (2.10) instead. 

3. Dark iV-soliton solution and its proof 

In this section, we show that the tau functions / and g representing the dark iV-soliton 
solution admit the compact determinantal expressions. This statement is proved by an 
elementary calculation using the basic formulas for determinants. We first prove that the 
proposed dark iV-soliton solution solves the bilinear equations (2.2) and (2.3) and then 
the trilinear equation (2.10) in place of (2.4). The implication of the equation (2.10) 
will be discussed in conjunction with the dark iV-soliton solution of the derivative NLS 
equation. Last, we perform the linear stability analysis of the plane wave solution (1.2) 
and provide a criterion for the stability. 
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3.1. Dark N -soliton solution 

The main result in this paper is given by the following theorem: 

Theorem 3.1. The dark N -soliton solution of the system of bilinear equations (2.2)-(2.4) 
is expressed by the following determinants 

f=\D\, (3.1a) 



9 = 



° * - 1 ° * (3.16) 



— 7* 1 



P 2 



z* t 



Here, D is an N x N matrix and z and z t are N -component row vectors defined below 
and the symbol T denotes the transpose: 

/ 7 \ , K - l Pi ( K P 2 , 1 > \ 

D = {d jk )i<j,k<N, d jk = b jk ^ —^ZjZ k , Zj=exp[pjX^ 1 -\ —^t + Q \, 

Pj +P k \ Pj Pj + IK J 

(3.2a) 

z = {z 1 ,z 2 ,...,z N ), zt=[ , J, (3.26) 

V Pi Pi Pn J 

where pj are complex parameters satisfying the constraints 

( Pj + w)(p* - in) = 1±^- Pj p*, j = 1,2,.., N, (3.2c) 

Kp 

Cjo (j — 1,2,..., N) are arbitrary complex parameters, 5j k is kronecker's delta and r is an 
auxiliary variable. 

The dark iV-soliton solution is parameterized by 2N complex parameters pj and 
Cjo (j — 1)2, ...,N). The parameters pj determine the amplitude and velocity of the 
solitons whereas the parameters (j determine the phase of the solitons. As opposed to 
the bright soliton case explored in I, however, the real and imaginary parts of pj are not 
independent because of the constraints (3.2c). Actually, it may be parameterized either 
by the velocity of the jth soliton or by a single angular variable, as will see in section 4. 
An auxiliary variable r introduced in (3.2a) will be used conveniently in performing the 
proof of (2.10). It can be set to zero after all the calculations have been completed. 

Remark 3.1. The tau function g given by (3.1b) is represented by the determinant of an 
(N + 1) x (N + 1) matrix. It can be rewritten by the determinant of an iV x N matrix. 
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To show this, we multiply the (N + l)th column of g by z* kt f p 
/cth column for k — 1, 2, N to obtain 



> 2 and subtract it from the 




(3.3) 



Although the tau function from (3.1b) is used in the proof of the dark iV-soliton solution, 
an equivalent form (3.3) will be employed in section 4 to analyze the structure of the 
solution. 

Remark 3.2. The complex parameters pj subjected to the constraints (3.2c) exist only 
if the condition k(1 + up 2 ) > is satisfied, as confirmed easily by putting pj = dj + ibj 
with real a,j and bj. We will show in section 3.6 that this condition is closely related to 
the stability of the plane wave solution of the FL equation. 

3.2. Notation and basic formulas for determinants 

Before entering into the proof of the dark iV-soliton solution, we first define the matrices 
associated with the dark iV-soliton solution and then provide some basic formulas for 
determinants. Although these formulas have been used extensively in I, we reproduce 
them for convenience. 

The following bordered matrices appear frequently in our analysis: 



where a, b, c and d are N component row vectors. Let Djk be the cofactor of the element 
djk- The following formulas are well known in the theory of determinants [10]: 




(3.4) 




(3.5) 



j,k=i 



D a T 
b z 



N 




(3.6) 



j,k=i 



D(a,b;c,d)||D| = |£>(a;c 



)|p(b;d)|-p(a;d)|p(b;c)|. 



(3.7) 
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The formula (3.5) is the differentiation rule of the determinant and (3.6) is the expansion 
formula for a bordered determinant with respect to the last row and last column. The 
formula (3.7) is Jacobi's identity. The proof of lemmas described below is based on the 
above three formulas as well as a few fundamental properties of determinants. 

3. 3. Differentiation rules and related formulas 

In terms of the notation (3.4), the tau functions / and g can be written as 

f=\D\, (3.8a) 

g=\D\ + ^\D(z* t ;z)\. (3.86) 
P 

The differentiation rules of the tau functions with respect to t and x are given by the 
following formulas: 

Lemma 3.1. 

/ t = i|L>(z t *;z)|-l|L>(z t *;z t )|, (3.9) 

/ x = -«|L>(z*;z)|+i|L>(z*; Z:c )|, (3.10) 
f xt = iK P 2 \D(z*; z)| - K \D(z* t ; z)| - k\D(z*; z t )| + i\D(z* t ; z x )\ 

-\D(z*, z* t ; z x , z)| + j\D[z\ z* t ; z, z t )\ - j 2 \D{z\ z* t ; z x , z t )\, (3.11) 

g t = i\D(z* t ;z)\ + \\D(z* tt ;z)l (3.12) 
P 

g x = i\D{z* t -z)\ + ^\D(z* t ;z x )\ + ^\D(z* t ,z*;z,z x )\. (3.13) 
Proof. We prove (3.9). Applying formula (3.5) to / given by (3.1) with (3.2a), one 
obtains 

N 

K — lpj 



ft = Kp > D jk —z j z i 

N I N 



P 

j,k=i j,k=i 
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where in passing to the second line, use has been made of the relation Zj ;t = (Kp 2 /pj)zj. 
Referring to formula (3.6) with z = and taking into account the notation (3.4), the 
above expression is equal to the right-hand side of (3.9). Formulas (3.10)-(3.13) can be 
proved in the same way if one uses (3.5), (3.6) and the relation z xt = np 2 z as well as some 
basic properties of determinants. □ 

The complex conjugate expressions of the tau functions / and g and their derivatives are 
expressed as follows: 

Lemma 3.2. 

/* = |L>|-i|L>(z*; z)|, (3.14) 
/; = -i|L>(z*;z t )| - ^\D(z*;z t )\ + -^|L>(z t *, z*; z t , z)|, (3.15) 

g* = \D\-i\D(z*;z)\ + ±\D(z*;z t )\. (3.16) 

Proof. It follows from (3.2a) that d* k = dkj + iz*Zk or in the matrix form, D* = D T + 
i(zjZl)f <jk<N . Since \D T \ = \D\, one has 



f* = \D + i( z j4h<j,k<N\ 



D z T 
-iz* 1 



Applying formula (3.6) to the right-hand side, formula (3.14) follows immediately. For- 
mulas (3.15) and (3.16) can be derived in the same way. □ 

3.4- Proof of the dark N-soliton solution 

3.4.1. Proof of (2.2) 
Let 

Pi = D t f-r-ip 2 (gg*-fr). (3.17) 

Substituting (3.8), (3.9) and (3.14)-(3.16) into (3.17), most terms are canceled, leaving 
the following three terms 

Pi = 4r{-|L>||L>(z t *,z*;z t ,z)| + |L>(z*;z)||L>(z t *;z t )| - |L>(z t *;z)||L>(z*;z t )|}. 
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This expression becomes zero by Jacobi's identity. □ 
3.4.2. Proof of (2.3) 

Instead of proving (2.3) directly, we differentiate (2.2) by x and add the resultant expres- 
sion to (2.3) and then prove the equation P 2 = 0, where 

P2 = f xt f* - Ut* - - /,/*) + kp 2 (99* - //*)■ (3.18) 

Substituting (3.8)-(3.11), (3.13) and (3.14)-(3.16) into (3.18) and rearranging, P 2 reduces 
to 

P 2 = ^{|D||D(z^z-z,zO|-|D(z-z)||D(z-z ^ )| + | J D(z*;z ^ )|| J D(z ^ *;z)|} 

+^{-\D\\D(z*,z; ] z x ,z t )\ + \D(z*;z x )\\D(z; ] z t )\-\D(z^ 

+\[-\D(z*-z)\\D(z\z;-z x ,z t )\+\D(z*-z x )\\D(zl 

9 (3.19) 
The first and second terms on the right-hand side of (3.19) vanish by virtue of Jacobi's 
identity. To show that the third term becomes zero as well, we consider the determinantal 
identity 



|L>(z*;z)| 


\D(z*;z)\ 


\D(z* t ;z)\ 




\D(z*;z x )\ | 


D(z*;z x )\ 


\D(z* t ;z x )\ 


= 


\D(z*;z t )\ 


D(z*;z t )\ 


\D(z* t ;z t )\ 





It is obvious that this determinant is zero since the first two columns coincide. The 
above assertion follows immediately by expanding the determinant with respect to the 
first column and using Jacobi's identity. Consequently, P2 = 0. □ 

Before proceeding to the proof of (2.10), we emphasis that the constraints (3.2c) have 
not been used in the process of the proof of (2.2) and (2.3). On the other hand, we find 
that the proof of (2.4) depends crucially on the constraints. This is an obstacle which has 
never been encountered in performing the proof of the bright 7V-soliton solution (see I). In 
conclusion, a direct proof of (2.4) still remains open and hence we shall prove the trilinear 
equation (2.10) instead. It turns out, however that its proof is found to be unfeasible. 
As we shall now demonstrate, introduction of an auxiliary variable r in the exponential 
function (3.2) would resolve this difficulty 
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3.4-3. Proof of (2.10) 



We first prepare the two lemmas to prove (2.10). The lemma 3.3 below gives a very simple 
relation between the partial derivatives f t and f T . It is to be noted that the constraints 
(3.2c) are used only for the proof of this lemma. 

Lemma 3.3. 



f t = (1 (3.20a) 

g t = (l + K P 2 )g T . (3.206) 

Proof. Extracting the factor Zj from the jih row and the factor z* k from the fcth column 

of the determinant \D\, respectively for j, k = 1, 2, N, one can rewrite the tau function 
/ into the form 



A? 



<j s jk + 



where 



K — ipj 

Pj + PkJ l<j,k<N 

Pj+P* 



Ci = (Pi + P*i) x + K P 2 [ — + ~Z )t + , 



T + Go + CjO> 



showing that / can be regarded as a function of Q (j = 1, 2, .., TV). Thus, differentiation 
of / with respect to t gives 



Similarly, one has 



N 



1\ df 



Kp f-M +p V9Q 



I * 

Pj + Pj 



df 



jr[ (Pj + ™)(p* -i/c) 90 ' 

The constraints (3.2c) are introduced into the above expression to give 

v 



/r 



+ 



n 9/ 



1 + «P 2 ^ VPi Pj* / 9 1 + ^ 



1 



This completes the proof of (3.20a). Repeating the similar procedure, one can show that 
the relation (3.20b) holds as well. □ 
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The lemma 3.4 gives the differentiation rules of / and g with respect to r: 
Lemma 3.4. 



/ T = i|£>(z»|, (3.21) 

/; = -i|£>(z*;z T )|, (3.22) 

9r = -^\D(z* t ;z)\ + hD(z* t ;z T )\, (3.23) 
K P P 

g XT = i|P>(z*;z)| + \\D(z* t ;z)\ + k\D(z*;z t )\ + -^\D(z* t ;z x )\ - ^|P>(z*; z r )| 
p z /tp z p z 

'P>(z*, z*; z, z x )| - 4l^K , z*; Zr, z)| + -^P(z t *, z*; z T , z x )|. (3.24) 



Kp z p z p z 



Proof. If one notes the relations 

i . 2 ^ . 

z tr — z t "I" IP z tj z it — z l/tZ T , 

which follows from the definition (3.2a) of Zj, the proof can be done straightforwardly 
along with the same procedure as that used in the proof of lemma 3.1 and lemma 3.2. □ 

With lemmas (3.2) and (3.3) at hand, we are now ready for starting the proof of (2.10). 

Proof of (2.10). If one replaces the t derivative by the r derivative in accordance with 
(3.20), the trilinear equation (2.10) can be rewritten in the form 

f*P 3 = f* P >, (3.25a) 
with the bilinear forms P 3 and Pg defined respectively by 

P3 = 9xrf ~ (fx ~ ™f)g T - -{g x f - gf x ), (3.256) 

K 

I'', //,/ ///, • ".•///• (3.25c) 

The trilinear equation (3.25) is proved as follows. Substituting (3.8), (3.10), (3.13), (3.23) 
and (3.24) into (3.25b) and applying Jacobi's identity to terms multiplied by \D\, P3 is 
simplified considerably. After some elementary calculations, one finds that 

P 3 = k|P(z*;z t )|{ \D\ + \\D(z* t ;z)\ - -^\D(z* t ;z x )\\. (3.26a) 
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Performing the similar calculation for P%, one obtains 

P 3 ' = i«{|£>| - \D(z*;z)\}{\D\ + ^\D(z* t -z)\ - ^\D(z* t ; z x )\ }. (3.266) 

Taking into account the formulas (3.14) and (3.22), the expressions (3.26a) and (3.26b) 
yield (3.25). The trilinear equation (3.25) coupled with lemma 3.3 now completes the 
proof of the trilinear equation (2.10). □ 

3. 5. Dark N-soliton solution of the derivative NLS equation 

In accordance with the fact that the FL equation is the first negative flow of the Lax 
hierarchy of the derivative NLS equation, the spatial part of the Lax pair associated 
with the former equation coincides with that of the latter equation with an identification 
q = u x [2,5]. This observation enables us to obtain the dark iV-soliton solution of the 
derivative NLS equation 

iqt + q xx + 2i(\q\ 2 q) x = 0, q = q(x,t), (3.27) 

under the boundary condition 

q pexp{i(KX - co't + ip {±) )} , x ±oo, (3.28) 

where u' = k 2 + 2np 2 and are real phase constants. In particular, we establish the 
following proposition: 

Proposition 3.1. The dark N-soliton solution of the derivative NLS equation (3.27) 
subjected to the boundary condition (3.28) is given in terms of the tau functions f and h 

by 

g = pe i(«*- w ' t )^ (3.29a) 

with 

f = \D\, h=\H\. (3.296) 
Here, D and H are N x N matrices defined respectively by 

D = {d jk )i<j,k<N, d jk = 5 jk ^ — ZjZ* k , zj = exp [pjX + {ip 2 - 2(k + p 2 )Pj}t + Qo] , 

Pj + Pk 

(3.30a) 
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H 



/ 7 \ 7 _ r K X Pj Pj * 
{lljk)l<j,k<N, I jk — Ojk - ~ — — — 2j2 fc , 



Pj + Pi Pi 

where pj are complex parameters satisfying the constraints 

p jP * = P 2 {k - i( Pj - p*)}, j = 1, 2, N, 
and (j (j = 1, 2, N) are arbitrary complex parameters. 



(3.306) 



(3.30c) 



Proof. The correspondence between q and u x mentioned above implies that the relation 



d 



pe 1 
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pe iKX j^(g x f -gfx + infg), 



dx V /, 

holds at t — 0. On the other hand, the expression in the parentheses on the right-hand 
side is just Pg defined by (3.25c) and hence it is equal to (3.26b). This fact and (3.14) 
lead, after applying the formula (3.6), to 

1 

p 2 ' " ' ' up" 

D 



g x f - gf x + wfg = m{\D\ - \D(z*; z)|}{|£>| + ^\D(z* t ; z)| - A^K; z x )|} 



i«r 



T 

i<i<JV 



Kk<N 



Multiplying the (N + l)th column of the determinant by KZ^/p* k and subtracting it from 
the fcth column for k = 1,2,...,N, one finds that the above expression becomes inf*h. 
Consequently, 



mpe 



f*h 



P 



t=o 



If one replaces q by iq and p by p/k, respectively and introduces the time dependence 
appropriately, one arrives at (3.29) with (3.30). The constraints (3.30c) follow from (3.2c) 
by the above replacement of p. The complex parameters pj subjected to the constraints 
(3.30c) exist only if the condition k + p 2 > is satisfied. □ 

It is instructive to perform the bilinearization of the derivative NLS equation under 
the boundary condition (3.28). This provides an alternative way to construct the dark 
iV-soliton solution given by proposition 3.1, as we shall see now. To this end, following 
the procedure used in [11, 12], we introduce the gauge transformation 



v cxp 



/X 
(P 2 
-oo 



\v \ 2 )dx 



(3.31a) 
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as well as the dependent variable transformation for v 

Then, equation (3.27) can be decoupled to the system of bilinear equations for / and h 

D x f-r-ip 2 (hh*-ff*) = 0, (3.32) 

D 2 J ■ f* - ip 2 D x h ■ h* + P \2k + p 2 )(hh* - //*) = 0, (3.33) 
\D t h ■ f + 2i(« + p 2 )D x h ■ f + D 2 x h ■ f = 0. (3.34) 
In view of (3.32), the modulus of v is given in terms of the tau function / by 

M 2 = p 2 + i|-hX, (3.35) 

which, combined with (3.31), yields the formula (3.29). Note from (3.31a) that \q\ 2 = \v\ 2 . 
It may be checked by direct computation that the tau functions / and h from (3.29b) 
with (3.30) satisfy the above bilinear equations. 

It is important to realize that we can take the limit k — > for the solution (3.29) since 
the dispersion relation is not singular at k — 0. This gives the iV-soliton solution of the 
derivative NLS equation on a constant background which has been studied extensively 
using various exact methods of solution such as the 1ST [13-16], Backlund transformation 
[17, 18] and Hirota's direct method [19]. On the other hand, for the dark iV-soliton 
solution given by (2.1), this limiting procedure is not relevant because of the singular 
nature of the dispersion relation. 

Last, we shall briefly describe the properties of the one-soliton solution for the purpose 
of comparison with those of the one-soliton solution of the FL equation. Introducing the 
new real parameters a\ and b\ by p\ = ai+i&i, the square of the modulus of the one-soliton 
solution from (3.29) and (3.30) with N — 1 can be written in the form 

i ,2 2 2a 2 sgna! 1 . . 

0i —p ; , , , N . 3.36a 



with 



0x = ai (x + Cl t) + 6 10 , ci = 2(6i + k + p 2 ), e 4<51 = a i + ^ + h ^ 



(3.366) 
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where sgnai denotes the sign of a±, i.e., a\ = 1 for ai > and a± = — 1 for a\ < 0, and 
9 W is a real constant. The constraint (3.30c) then becomes 



a\ + b{ = p 2 (2b 1 + K). (3.37) 

Using (3.36b) and (3.37), the parameters a± and b\ are expressed in terms of the velocity 
Ci of the soliton as 

a\ = ^ (c m ax - ci) (ci - c min ) , 61 = ^ - « - p 2 , c min < Ci < c max , (3.38a) 
where 

c m ax = 2(k + 2p 2 ) + 2p v / /« + P 2 , c min = 2(k + 2p 2 ) - 2py / K + p 2 . (3.386) 

One must impose the condition k + p 2 > to assure the existence of soliton solutions. 
Recall that this condition coincides with a criterion for the stability of the plane wave 
(3.28) [20]. We see from (3.36) that if ai > 0, then \qi\ takes the form of a dark soliton 
whereas if a\ < 0, it becomes a bright soliton on a constant background u = p. 

Let Ad and At, be the amplitudes of the dark and bright solitons, respectively with 
respective to the background. The amplitude-velocity relations follow from (3.36) and 
(3.38). They read 

A d = p- y/ Cl - k - 2p 2 - + p 2 , (3.39a) 

A b = xAi - k - 2p 2 + y/ K + p 2 - p. (3.396) 

The detailed analysis for the case k > has been undertaken in [21]. To sum up, the 
solution has been shown to exhibit the spiky modulation of the amplitude and phase. 
It also has been demonstrated that the bright soliton reduces to an algebraic soliton for 
both limits c\ — > c max and c\ — > c mm whereas the algebraic dark soliton never exists. In 
the case k < which has not been treated in [21], however, a careful inspection of (3.36) 
and (3.38) reveals that the algebraic bright and dark solitons are produced in the limit 
ci — > c max and c\ — > c m i n , respectively. The latter new feature is pointed out here for the 
first time. 
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Remark 3.3. Using the result obtained in proposition 3.1, we can construct the dark 
iV-soliton solution of the modified NLS equation 



iqt + Qxx + v\q\ 2 q + i7(l<?| 2 <?):r = 0, q = q(x, t), 



(3.40) 



under the boundary condition 

q — > pexp {i [kx — u"t + } , x — > ±oo, 



(3.41) 



where oj" = k 2 — pp 2 + •ynp 2 and p and 7 are real constants. To show this, we apply the 
gauge transformation 

2 



q = exp 



7 



7 



_ 2/i - 
g, x = x H t, t — t, 

7 



(3.42) 



to equation (3.40) and see that it can be recast to the derivative NLS equation iqi + q xx + 
7(l<?| 2 <?):r = 0, which coincides with equation (3.27) with the identification q = q, x — 
x, t — t and 7 = 2. The dark iV-soliton solution of the equation (3.40) then takes the 
form 

q = pe i ^" t) ^, (3.43a) 
where the tau functions /' and h! are given respectively by 



E 

7 



K — — — lpj 



ti 



Pj + Pi 



Pl+Pl Pi '" 



Z j Z k 



l<j,k<N 



l<j,k<N 



with 



(3.436) 



(3.43c) 



(3.43d) 



Zj = exp [pjx+{ip 2 - (2k + -fp 2 )pj}t + Co] • 
The constraints for pj become 

PiPi = ^{«-^-i(Pi-Pi)}, j = l,2,...,N. (3.44) 
The complex parameters pj exist only if the condition 7 — ^ + ^ j > is satisfied. 
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The following two special cases are worth remarking. The case p = and 7 = 2 
reduces to the result given by proposition 3.1. On the other hand, in the limit 7 — > 
while p being fixed, we first replace Zj by y/^jZj for j = 1,2, ...,N and then take the limit, 
producing the dark iV-soliton solution of the NLS equation. Note, in this limit, that the 
constraints (3.44) reduce to pjp* = — pp 2 /2 and hence the dark soliton solutions exist only 
if the condition p < is satisfied. 

3. 6. Stability of the plane wave 

We have considered the dark solitons on the background of a plane wave pe 1 '" 1 ""'' with 
ijj = 1/ ' n + 2p 2 . It is important to see whether the background field is stable or not against 
perturbations. If unstable, then dark solitons would not exist, as will be demonstrated in 
the next section. To this end, we perform the linear stability analysis of the plane wave. 
Following the standard procedure, we seek a solution of the form 

u = (p + Ap) e i ( ra -^+ A *) , (3.45) 

where Ap = Ap(x,t) and A<p = A(f)(x,t) are small perturbations. Substituting (3.45) 
into the FL equation (1.1) and linearizing about the plane wave, we obtain the system of 
linear PDEs for Ap and A<p 

Ap xt + p(u - 2p 2 )A(p x - npA(j) t - Anp 2 Ap = 0, (3.46a) 

pA<t> xt - (w - 2p 2 )A Px + kA Pi = 0. (3.466) 

Assume the perturbations of the form e 1 ^-^) with A real and v possibly complex and 
substitute them into (3.46) to obtain a homogeneous linear system for Ap and A0 

(\v - 4/tp 2 ) Ap + i{p\(u - 2p 2 ) + /tpz/}A0 = 0, (3.47a) 

-i{(w - 2p 2 )A + kv}A p + pAz/A0 = 0. (3.476) 
The nontrivial solution exists if v satisfies the quadratic equation 

A 2 

(A 2 - k 2 )u 2 - 2(2np 2 + l)\u - — = Q. (3.48) 
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Solving this equation, we obtain 



A 

v 



\ 2 -K 2 

Thus, if the condition 



2np 2 + 1 ± -^\ 2 + Ak 3 (k P 2 + l)p 2 

K 



(3.49) 



k(k P 2 + 1) > 0, (3.50) 



is satisfied, then v becomes real for all values of real A, implying that the plane wave is 
neutrally stable. It is evident that this condition always holds for k > 0. For negative 
k, on the other hand, we put k = —K with K > and see that the stability criterion 
turns out to be as Kp 2 > 1. Last, we remark that a similar stability analysis has been 
performed recently in conjunction with a plane wave solution of the original version of 
the FL equation [22, 23]. 

4. Properties of the soliton solutions 

In this section, we detail the properties of the soliton solutions. To this end, we first 
parametrize the complex parameters pj and Q by the real quantities aj, bj, 9jo and Xjo as 

Pj = a j + Go = djo + iXjo, 3 = !» 2 > •••> N , (4- 1 ) 

and introduce the new independent variables 9j and Xj according to the relations 

0j = dj(x + Cjt) + 6j , cj = 2 jo , j = 1, 2, N. (4.2a) 

Xj = bj(x - Cj t) + X jo, J = 1,2,..., N. (4.26) 
In terms of these variables, the variables Zj defined by (3.2a) are put into the form 



z 



= e 9 ^, j = l,2,...,N, (4.2c) 



after setting r = 0. Substituting (4.1) into (3.2c), the constraints for pj can be rewritten 
as a quadratic equation for bj 

b 2 - 2K 2 p% + a 2 - K 3 p 2 = 0, j = 1, 2, .., N. (4.3) 
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The solution to this equation is found to be as follows: 

bj = M 2 ± W) -4 J = l,2,...,N. (4.4) 

We can see from the above expression that the real bj (j = 1,2, ...,iV) exist only when 
the condition /t 3 p 2 (l + /tp 2 ) > is satisfied. This coincides with the criterion (3.50) for 
the stability of the plane wave, as discussed in section 3.6. Throughout the analysis, we 
assume this condition to assure the existence of soliton solutions. It is to be noted from 
(4.2) and (4.3) that the parameters a,j and bj are expressed in terms of Cj as 

K 2 , _ , . 1 



2 



(Cmax Cj) (Cj C m in) i bj 2/CC ■ ^ ^ ^j)' C m ; n <! Cj <! C max , (4.5d) 



where 

c m ax = ^ {l + + 2 v //tp 2 (l + / t p 2 )} , c min = {l + 2/tp 2 - 2^/^(1 + /tp 2 )} . 

(4.56) 

The relations (4.5) correspond to (3.38) for those of the one-soliton solution of the deriva- 
tive NLS equation. Thus, the dark iV-soliton solution is characterized by the N velocities 
Cj {j — 1, 2, N) and the 2N real phase constants 9j and Xjo (j = 1, 2, JV), the total 
number of which is 3./V". 

Another parameterization of the solution is possible if one introduces the angular 
variables 7j - by 

a, = a/ /t 3 p 2 (l + /tp 2 ) sin 7j, (4.6a) 
bj = (k P ) 2 + v^y(i + Kp 2) cos 7j> o < 7j < 2vr, 7j ^ tt, j = 1, 2, TV. (4.66) 
In terms of 7j -, p^ from (4.1) can be written in the form 

p, = i{M 2 + v^VOT^ye^}, j = l,2,...,iV, (4.7) 
and the velocity Cj of the jih soliton given in (4.2a) is expressed as 

Cj = / . 2 x I 1 + 2k P 2 -2sgn/tv//tp 2 (l + /tp 2 ) cos 7j | . (4.8) 

It follows from the above parametric representation that pj lies on the circle of radius 
a/ /t 3 p 2 (l + up 2 ) centered at i(/tp) 2 in the complex plane. Let us first describe the proper- 
ties of the one- and two-soliton solutions and then address the general iV-soliton solution. 
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4-1- One-soliton solution 

The tau functions / = f\ and g — g\ for the one-soliton solution follows from (3.1)-(3.3) 
with N = 1. They read 

/i = l + : — 0i = 1 : — *~ 4.9 

Pi +Pl Pi+ Pi Pi 

The one-soliton solution ui follows from (2.1) with (4.9), yielding 

-y _ K+6i+iai ai+ibi g 20i 

Ul = pe^^-^ 2a ! ; ai -' l \ . (4.10) 

The above expression can be put into the form 

Ul = \ Ul \ e i{KX ~ ut) exp {i ((f) + (+) ) } , (4.11) 

where the square of the modulus of U\ is represented by 
i ,2 2 2a?csgn(/tai) 1 



y/al + + h) 2 cosh 2(0i + 5i) 



(jt+bi) sgnai 
v / a 2 + ( K+fel )2 



with 



0i = "i^ + '-iO + 9 W . Cl = e 4 * = a ' + ( ; 2 +&l)2 , (4.126) 



and the tangent of the phase and (j)^ being given respectively by 

^_ {a? + &!(« + &!)} cosh2(g 1 + ^ 1 ) + 5 1 sgna lv /^ + (/t + & 1 ) 2 
tan0_ «a! sinh 2^ + 5,) ' 



(4.136) 



tan 0« = + + M 

It can be confirmed by direct substitution that (4.11) indeed satisfies the FL equation. The 
one-soliton solution (4.10) is a one-parameter family of solutions. The parameterization 
in terms of a\ will be employed in classifying the soliton solutions. The parameters C\ 
and b\ are then expressed by a\. See (4.2a) and (4.4) whereas the parameters p and k are 
fixed by the boundary condition (1.2). The relation (4.5) will be used conveniently when 
considering the generation of algebraic solitons in the limit |ai| — > 0. The form of \u±\ 
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from (4.12) reveals that If na x > 0, then \ui\ takes the form of a dark soliton whereas if 
Kdi < 0, it becomes a bright soliton on a constant background u = p. Note from (4.12) 
that the width of the soliton may be defined by (2|ai|) _1 . The net change of the phase 
caused by the effect of nonlinear modulation is given by (4.13). Roughly speaking, the 
phase behaves like a step function as a function of 9\. Specifically, a rapid change of 
the phase occurs in the vicinity of the center position of the soliton {6\ = —5i), yielding 
a phase difference % (or — n). As a result, the phase of U\ changes by a quantity 20( + ) as 
9i varies from — oo to +oo, where (p^ is given by (4.13b). 

Let us classify the one-soliton solutions in accordance with the sign of k. We consider 
the two cases, i.e., case 1 (k > 0, ai «s 0) and case 2 (n < 0, ai ^ 0) separately. For each 
sign of k, both dark and bright solitons arise, as we shall show now. 

4.1.1. Case 1: k > 

In this case, the velocity c\ of the soliton is positive, as evidenced from (4.12b). Let A d 
and A b be the amplitudes of the dark and bright solitons, respectively with respective to 
the background. We then find from (4.5) and (4.12) that 



A d = p 



y / p 2 -2c 1 { V / "a? + ( K + 6 1 ) 2 -(^ + &i)} 



1 



y/c - a/1 + np 2 , ai > 0, ci = c> 0, (4.14) 



A b = ^p 2 + 2 Cl \^Ja\ + (k + 6i)2 + (k + 6i)J - p 

= (nyfc + a/1 + np 2 ] - p, ai < 0, Ci = c> 0, (4.15) 

where c = |ci| lies in the interval c m i n < c < c max with c max and c m i n being given by (4.5b). 
Note from (4.5a) that k + b\ = (1 + k 2 c\)/ (2kci) > for k > and ci > 0. This estimate 
will be used to judge the existence of algebraic solitons in the limit of infinite width. 

Figure 1 plots the dependence of the amplitudes A = A d and A = A b on the velocity 
c = |ci| for p — 1 and k — 2. 

(i) Dark soliton: a x > 
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< 1.5 




Figure 1. Amplitude- velocity relation for the dark soliton Ad (solid line) and bright 
soliton Af, (broken line) for p — 1 and k — 2. 



As seen from figure 1, the amplitude Ad of the dark soliton becomes an increasing function 
of the velocity c in the interval c m i n < c < cq and a decreasing function in the interval 
c < c < c max , where c max (71 = ir) and c min (7! = 0) are given by (4.5b) and a critical 
velocity Cq and the corresponding angle 70 by 



c = 



1 + up 1 



K 



at 71 = 70 = cos 1 



(k P 2 )^(3 + 2k P 2 ) 



(0<7o<tt)- (4.16) 



2{l + K P 2 )i 

In the present numerical example (p — 1,k — 2), c m j n = 0.025, Co = 0.75, c max = 2.47. 
The above observation shows that in the interval c < c < c max , a small dark soliton 
propagates faster than a large dark soliton. A similar behavior has also been found in I 
for the bright soliton solutions of the FL equation with zero background. 

Figure 2 depicts the profile of U — \ui\ at t — for three different values of c, i.e., a: 
c = c = 0.75(7i = 70 = 0.90tt), b: c = 0.33(7i = 5tt/6), c: c = 0.098(7i = 2vr/3) with 
the parameters p — 1, k — 2, 6 W — —Si and xw — 0. When c = c , the amplitude of the 
dark soliton attains the maximum value Ad = p. See figure 2 a. It then turns out that 
the intensity of the soliton center falls to zero. Such a soliton is well-known in the field 
of nonlinear optics. It is sometimes called a black soliton. For this specific value of c, one 
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Figure 2. Profile of the amplitude of the dark soliton U — \ui\ at t — 0. a: c = Cq = 0.75, 
b: c = 0.33, c: c = 0.098. The profile a is a black soliton. 



finds from (4.5), (4.12) and (4.13) that 



ai = 



«2p(4 + 3/tp 2 )2 



M 2 



(4.17a) 



2(1 + /tp 2 ) 



2(l + /tp 2 )' 





tan <p = tanh(^! + Si), 



tan 



6 



(4.176) 



The profile of |wi| 2 



2 from (4.12) then becomes 



|2 2 i 



4 + 3/tp 2 1 



(4.18) 



2(1 + up 2 ) cosh 2(0i + + 




As confirmed easily from the above expression, the minimum value of |wi| is zero at 
Q\ = —5\. The algebraic dark soliton may be produced from (4.12) by taking the limit 
ai — > +0. However, as already noticed, the value of k + bi is positive so that \ui\ tends 
simply to a constant value p. Hence, this limiting procedure is irrelevant for the dark 
soliton solution under consideration, indicating that the algebraic dark soliton does not 
exist for k > and a\ > 0. 

Figure 3 shows the profile of Mr = Re[wi] at t — 1 for the black soliton. The broken 
line indicates ±|«i| (see figure 2 a). One can see that the dark soliton exhibits phase 
modulations near the center position of the soliton. This peculiar feature is in striking 
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Figure 3. Profile of a black soliton Mr = Reu± at t — 1. 



contrast to the bright soliton solution of the NLS equation for which no phase modulation 
occurs. A similar behavior has been observed for both dark and bright soliton solutions 
of the derivative NLS equation with the background of a plane wave [21, 24]. 

(ii) Bright soliton: a x < 

Figure 4 depicts the profile of the bright soliton U — \ui\ at t — for three different values 
of c, i.e., a: c = 2A7(^ = l.OObr), b: c = 0.73(7i = l.br), c: c = 0.025(7i = 1.999tt) 
with p = 1 and k — 2. The feature of the bright soliton differs substantially from that 
of the dark soliton. To be specific, the amplitude of the bright soliton always becomes 
an increasing function of the velocity (see figure 1). It takes the maximum value at 
c = c max (7i — > 7T+0, ai — > —0) and the minimum value at c = c min (7i — > 2n — 0, a\ — > — 0). 
At these limiting values of the velocity, the algebraic soliton is produced from the soliton 
of hyperbolic type. Indeed, if we put 9 W = aix — 5i in (4.10) and (4.12) with xo being a 
real constant and then take the limit ai — > —0, we find 

Ul = pe'(— ) X + Ct + Xo ~ i ^S k : (4 . 19a) 
x + ct + x - i 2{k+6i) 

| Ul |2 = p 2 + 1 {4m) 

l + *»c (x + ct + Xo)a+ (_«,_)»' 
where 6i = (1 — k 2 c) /2kc by (4.5a) and c = c max or c min . Note from (4.12b) that 6^ = K,p 2 /c 
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Figure 4. Profile of the amplitude of the bright soliton U — \ui\ at t 
b: c = 0.73, c: c = 0.025. The profiles a and c are algebraic solitons. 



0. a: c = 2.47, 



when ai — > —0. One can see that the algebraic soliton has no free parameters except a 
phase constant xq since the velocity c is determined by p and k which are fixed by the 
boundary condition. 

To derive (4.19a) from (4.10), we use the following expansion formulas for small a\\ 

2 l a l| j ai (x+ct+x ) 2 l a ll 



v / a 2 + (/t + fei)2 

k + bi — iai 



-{l + 2a 1 (x + ct + x ) + <3(a 2 )}, (4.20a) 



2ai 



e 26 * 1 ~ sgnai sgn(/t + 6i) 



Ik + 6 



1 + ai < 2(x + ct + xq) — i 



1 



« + h 



+ 0(aj) 



k + 6i + iai ai + i6i 



,20i 



2ai 



ai — i6i 



-sgnai sgn(/t + 6x) 



1 + a-i < 2(x + ct + x ) — i 



2/t + 6i 



+ 0(a 2 ) 



(4.206) 



(4.20c) 



6i(k + 6i) 

Because of the inequalities a\ < and « + bi > in the current problem, one finds 
that the condition sgn ai sgn(/t + & x ) = —1 is satisfied, which yields (4.19a) by taking the 
limit ai — > —0 for (4.10). Actually, under the above condition, the leading-order terms of 
the denominator and numerator of (4.10) turn out to be of order a±. Consequently, the 
expression (4.10) has a limiting form (4.19a) in the zero limit of a±. On the other hand, 
the expression (4.19b) follows either directly from (4.19a) or from (4.12) by performing 
the similar limiting procedure. 
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Figure 5. Profile of an algebraic bright soliton ur = Re^i at t — 1. 



A representative profile of the algebraic bright soliton U = \ui\ at t — and the 
corresponding profile of % = Re«i at i = 1 are shown in figure 4 a and figure 5, 
respectively. 

The novel feature of the bright soliton mentioned above deserves a few comments. 
First, the amplitude of the bright soliton tends to a finite value when its width tends to 
infinity, as opposed to the behavior of the dark soliton discussed just before for which the 
amplitude becomes zero in this limit. Second, the FL equation has an infinite number of 
conservation laws [3]. Among them, we evaluate the conserved quantity I = f^duxl 2 — 
n 2 p 2 )dx for the one-soliton solution (4.10). This quantity may be termed the energy of the 
soliton in accordance with the correspondence between the solution u of the FL equation 
and the solution q of the derivative NLS equation. Using the relation (|M x | 2 )t = (M 2 ) x 
which follows directly from the FL equation, we obtain 

I = —4 sgn ai tan 

We find from this expression that in the limit of infinite width |ai| — > 0, I becomes zero 
for the dark soliton (a x > 0) and tends to a finite value 2ir for the bright soliton (ai < 0). 
See also an analogous calculation for the bright soliton solution of the derivative NLS 
equation with zero background [25]. 

4.1.2. Case 2: k < 



1 



a\ + (k + h) 2 — (k + &i)sgnai 
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Figure 6. Amplitude-velocity relation for the dark soliton A d (solid line) and bright 
soliton At, (broken line) for p = 1 and k = —2. 



For negative k, the expressions of the amplitude for the dark and bright solitons are given 
respectively by 



A d = p-^jp 2 + 2c 1 |y / af + (/t + 6 1 ) 2 + (k + &i)J 

K^~c - ^Kp 2 - 1 , ai < 0, ci = -c < 0, (4.21) 



= P 



A b =^jp 2 -2c i yaj + (K + h) 2 - (k + &0 j - p 

= -L (V^ + V^P 2 - l) - P, ai > 0, ci = -c < 0, (4.22) 
y K ^ ' 

where K = —k is a positive wavenumber and the velocity c lies in the interval c' min < c < 

C with 

4ax = ^ {2^P 2 -1+ 2 V ^W Z 1)} , 4in = ^ {2^P 2 -1-2v/^W^T)} . 

(4.23) 

Recall that the condition Kp 2 — 1 > must be imposed to assure the existence of the 
soliton solutions. 

Figure 6 plots the dependence of the amplitudes A = A d and A = A b on the velocity 
c = | Ci | for p = 1 and k — —2. When compared with figure 1 for k > 0, there appear 
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Figure 7. Profile of the amplitude of the dark soliton U — \u\\ at t — 0. a: c = Co = 0.25, 
b: c = 0.16, c: c = 0.043. The profile a is a black soliton and the profile c is an algebraic 
soliton. 

several different features for k < 0. In particular, the algebraic dark soliton would arise 
in the limit c — > c' min since in this limit, the amplitude A d tends to a finite value. In 
addition, the algebraic bright soliton exists only in the limit c — > c' m3uX . We now proceed 
to the detailed description of the soliton solutions. 

(i) Dark soliton: a x < 

It follows from (4.5) with k = —K, c\ = —c that k + b± — \/2Kc — K/2. Since c' min < c < 
c max by (4.23), the possible value of k + b\ is restricted by the inequality 

< k + bi < K 



K 



K p 



y/Kf?{K? - 1) 



Kp 2 - 1 + ^Kp 2 {Kp 2 - 1) . (4.24) 



One can see that the upper limit of k + b\ is attained when c = c' min and its limiting value 
is positive by the condition Kp 2 > 1 whereas the lower limit is attained when c = c' max 
and is negative. In view of this fact, the algebraic dark soliton would be produced in 
the limit c — > c' min for which sgn(/t + bi) > 0. Actually, taking the limit a\ — > —0 for 
the solutions (4.10) and (4.12) and using the expansion formulas (4.20), we find that the 
hyperbolic soliton reduces to the limiting form 

x cl -t- XQ 1 2b 1 (-K+b 1 ) 



u\ = pe 



i(-Kx-ut) 



x ct + x i 2 (_^ +bi ) 



(4.25a) 
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Figure 8. Profile of an algebraic dark soliton u R = Re Ui at t — 1. 



Ml 



P 2 - 



2Kc 2 



,2' 



(4.256) 



1-K*c ( x - ct + x ) 2 + (^j) 1 

where &i = (1 + K 2 c)/2Kc and c = c^ in . Since 1 — Ifc^ > by virtue of the condition 
Kp 2 > 1, the expression (4.25b) actually represents an algebraic dark soliton. 

The black soliton appears when the velocity c takes a specific value c = c' , where 



c' Q = (Kp 2 - 1)/K 2 at 7i = i Q = cos' 1 



(Kp 2 )H3-2Kp 2 ] 



2(Kp 2 - 1) 



Its profile is represented by 



Pi 



P 



3Kp 2 



(n<i < 2tt). 

(4.26) 

(4.27) 



2(^P 2 -l)cosh 2(9, + 5,) + ^^ 

It is important to notice that the inequality c' min < c' Q < c' max requires the condition 
Kp 2 > 4/3 for the wavenumber K. It then turns out that expression (4.27) takes the 
form of a black soliton. 

Figure 7 depicts the profile of U — \ui\ at t — for three different values of c, i.e., a: 
c = c' = 0.25(7i = 7o = 5vr/4), b: c = 0. 16(71 = Air/ 3), c: c = 0.043(7i = 2tt) with the 
parameters p — 1, k — —2, 9 W = —Si and xio — 0. In this example, c' min = 0.043, c' Q = 0.25 



and cL 



1.46 (see figure 6). An algebraic soliton appears at the lower limit of the 
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x 

Figure 9. Profile of the amplitude of the bright soliton U — \ui\ at t — 0. a: c = 1.46, 
b: c = 0.81, c: c = 0.19. The profiles a is an algebraic soliton. 

velocity, i.e., c = c' min whereas a black soliton arises at c = c' . Figure 8 shows the profile 
of ur = Re u\ at t = 1 for an algebraic dark soliton. 

(ii) Bright soliton: a\ > 

The crucial difference between the case 1 and the case 2 for the bright solitons is observed 
if one compares figure 6 with figure 1. Notably, the bright soliton with k < reduces to 
an algebraic soliton only at the upper limit of the velocity c = c' max whereas the bright 
soliton with k > has two critical velocities c max and c m i n for which algebraic solitons 
are produced. Figure 9 depicts the profile of U — \u±\ at t — for three different values 
of c, i.e., a: c = 1.46(71 = 0.998tt), b: c = 0.73(7i = 0.9tt), c: c = 0.025(7i = 0.7tt) with 
p = 1 and /t = —2. Figure 10 shows the profile u R = Reui of an algebraic bright soliton 
at t — 1 which corresponds to the profile a in figure 9. 

4-1.3. Note on algebraic solitons 

We have seen that the algebraic solitons arise from the hyperbolic solitons when certain 
conditions are satisfied. Here, we summarize the result. The algebraic bright solitons 
are produced when the conditions sgnaisgn(ft + bi) = —1 and sgn(«;ai) < are satisfied 
simultaneously whereas the corresponding conditions for the dark algebraic solitons are 
given by sgnaisgn(/t + bi) = —1 and sgn(Kai) > 0. Thus, for k > 0, the conditions 
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Figure 10. Profile of an algebraic bright soliton ur = Reui at t — 1. 

sgn(/t + 61) = 1 and sgn(/t + bi) = —1 are responsible for the generation of the algebraic 
bright and dark solitons, respectively. Since k + b\ > in this case, only the bright 
algebraic soliton exists. See figure 1. For k < 0, on the other hand, the above conditions 
turn out to be sgn(/t + bi) = —1 and sgn(/t + 61) = 1, respectively. Under this setting, 
the limiting value of k + b\ becomes negative for the bright soliton and positive for the 
dark soliton, respectively, implying the existence of both types of algebraic solitons. See 
figure 6. In conclusion, we emphasize that the criterion for the existence of solitons (which 
depends crucially on the sign of k) plays an important role in our analysis. 

4-2. Two-soliton solution 

As clarified by the analysis of the ono-soliton solutions, both dark and bright solitons 
exist in our system. Therefore, the two-soliton solutions can be classified into three types, 
i.e., dark-dark solitons, dark-bright solitons and bright-bright solitons. Here, we focus our 
attention on the dark-dark solitons. Especially, we investigate the asymptotic behavior of 
the solution for large time. The two-soliton solution describing the interaction between a 
dark soliton and a bright soliton will be briefly discussed. For both cases, we assume that 
k > 0. 

4-2.1. Dark- dark solitons 

The tau functions f 2 and g 2 representing the dark two-soliton solution are given by (3.1)- 
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(3.3) with N = 2 subjected to the conditions re > 0, a x > 0, a 2 > 0. They read 

ipi)(re-ip 2 )(pi -P2)(p{ -Pt 



7 Z 2 Z 2 



* i 



(Pi + Pi) (Pi + P%) (P2 + P*) (P2 + P2) 

(4.28a) 

(k + ip\){K + ip^)(pi - P2)(pj - P2) P1P2 

Z\Z 2 Z-yZ 2 . 



h = H ; — j + - 

Pi + P* P2 + P 2 

_ K + igi Pl , « + IP2 P2 
^ Pi + P\ P\ P2 + P2 P2 ^ ' (Pi + Pi ) (Pi + P2) (P2 + Pi) (P2 + P2) PlP2 

(4.286) 

To investigate the interaction process of two solitons, we first order the magnitude of 
the velocity of each soliton in the (x,t) coordinate system as c\ > c 2 > 0. Invoking the 
definition (4.2a) of the velocity of the solitons, this can be established by imposing the 
condition \p ± \ < \p 2 \ on the amplitude parameters. Now, we take the limit t — > — 00 with 
61 being fixed. Since in this limit \z x \ =finite and \z 2 \ — > 00, the leading-order asymptotics 



of f 2 and #2 are found to be as 
re - ip 2 



P. 



92 



P2 + P2 

re + ip 2 p 2 



z 2 z* 2 < 1 + 



P2 + P2 P2 



^2^2 



(/c-ipi)(pi -p 2 )(Pi -P2) 

(pi +P* 1 )(pi+P* 2 )(P2+P*l)) 

(K + ip*)(Pi -P2W1 



p 2 ) Pi * 

7 ^1 



(4.29a) 



(4.296) 



(Pi +P*i)(Pi +P* 2 )(P2 + P*)p* 
The asymptotic form of the two-dark soliton solution follows from (2.1) upon substituting 
(4.29) into it, giving rise to 



u 2 ~ p exp ji {^Kx — ut + <f>[ ^ I — 



Pi ^/ z t * 
Pi+Pj Pi 1 1 

h Pl+P* 



2-i 2 



1*1 



where 



z\ exp 



X ^ = ar£ 



-In 



Pi +P2 



Pi -P2 

re + ips P2 



+ 7T. 



(4.30a) 



(4.306) 



(4.30c) 



- J P2 P2, 

Let Mi(^i) be the dark one-soliton solution (4.10). Then, the asymptotic form of u 2 can 
be written in terms of Mi as 



U2 



exp (i^ _) ) Ml (^ + A^ } ), A^ _) = -In 



Pi +P2 



Pi -P2 



(4.31) 
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Next, we take the limit t — > +00 with Q\ being fixed. In this limit, \z\\ =finite and 
\z2 1 —7- 0. Therefore, the tau functions f 2 and g 2 and the two-soliton solution w 2 behave 
like 

(4.32) 



Pi + p* 



K + ip x pi 



Pi + Pi Pi 



«2 ~ pe 



l(KX-Ult) 



1 _ K+ipi pi * 
1 Pl+pj Pi 

P1+P1 1 1 



It follows from (4.33) that 



u 2 ~ Ul (0! + A^ +) ), A#f ; = 0. 



3 (+) 



(4.33) 



(4.34) 



The trajectory of the center position x = x c (t) of the jth soliton is described by the 
equation 0j + A0^ = 0, or x c = —Cjt — (9j + AQ^)/a,j. Since the soliton propagates to 
the left, the phase shift Axj of the jth soliton can be defined by the relation 



Axj = x c (-oo) - x c (+oo) = i- (A0<. +) - A0<- _) ) , j = 1, 2. 
We see from (4.31) and (4.34) that the fast soliton suffers a phase shift 



(4.35) 



Axi = — In 

Oi 



P1-P2 



(4.36) 



In terms of the angular variable 71 and 72 defined by (4.6) and (4.7), this expression can 
be rewritten in the form 

sin \ (71 + 72) 



Ax\ = — In 

Oi 



sin \ (71 - 72) 



y/K 3 p 2 (l + up 2 ) sin 71, < 71 < vr. (4.37) 



We can perform the similar asymptotic analysis while keeping 62 fixed. Hence, we 
quote only the final results. As t — > —00, the expressions corresponding to (4.29) and 
(4.31) read respectively 



f ,,K-\p 2 
h ~ 1 + —7 Z 1 Z 2 



K + \p* 2 p 2 



P2 + P2 P2 + P2P2 

u 2 ^ Ul (0 2 + A6 { f ] ) , A0£ _) = 0. 



(4.38) 
(4.39) 
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Figure 11. The interaction of two dark solitons. 



As t — > oo, on the other hand, they take the form 
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Pl + Pl 

k + ipi pi 



z x z\ \ 1 + 



(/C-jp 2 )(Pl -P2)(Pi~P2) 
(P2+P^)(Pi+P^)(P2+Pi)) 



Z 2 Z 2 



Pi + Pi Pi 



^ I (P2+P 2 )(Pl+P 2 )(p 2 +Pl))p 2 ^ 

m 2 ~exp (i0< +) ) Ul (0 2 + A0< +) ), 



AO. 



(+) 



-In 



P2+P1 



P2 ~P\ 



4 +) = arg 



k + ip* pi 
« - ipi Pt 



+ 7T. 



(4.40a) 

(4.406) 
(4.41a) 
(4.416) 



The phase shift of the slow soliton follows from (4.35), (4.39) and (4.41), resulting in 



Ax 2 = In 

a 2 



P2+PI 



P2-P1 



(4.42) 



or equivalently in terms of the angular variables 71 and 72, it reads 



Ax 2 = In 

a 2 



sin \ (72 + 71) 



sin \ (72 - 71) 



a 2 



= a/k 3 p 2 (1 + np 2 ) sin 72, < 72 < n. (4.43) 



An inspection of the formulas (4.36) and (4.42) reveals that Axi > and Ax 2 < under 
the setting a\ > 0, a 2 > 0. 
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Figure 12. The interaction between a dark soliton and a bright soliton. 

Figure 11 shows the intercaction of two dark solitons with the parameters p = 1, k — 
2, a = 0.75(7i = 0.907r),c 2 = 0.24( 72 = 0.80tt) and Cio = C20 = so that from (4.14), 
A dl = 1.0 and A d2 = 0.47. It can be seen from figure 1 that the amplitude of each dark 
soliton is an increasing function of the velocity for the present choice of the parameters. 
Note, in this example, that the large soliton is a black soliton since its asymptotic ampli- 
tude is A d i — P — 1- The phase shifts evaluated from the formulas (4.37) and (4.43) are 
given by Axi = 0.70 and Ax 2 = —0.36, respectively. Figure 11 shows clearly a typical 
interaction process of solitons, i.e., as time goes, the large soliton gets close to the small 
soliton and overtakes it and after the collision, both solitons eventually separate each 
other without changing their profiles. The net effect of the collision is only the phase 
shift. 

4-2.2. Dark-bright solitons 

The two-soliton solution consisting of a dark soliton and a bright soliton is obtained by 
choosing the parameters such as k > 0, a x > and a 2 < 0, for example. The asymptotic 
analysis can be performed as well for this solution and hence the detail will be omitted. 

Figure 12 depicts the interaction between a dark soliton and a bright soliton with 
the parameters p = = 2, ci = 0.75(7i = 0.907r),c 2 = 0.24(7 2 = 1.2n) and Cio = 
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C20 = 0, showing that the dark soliton propagates faster than the bright soliton. The 
asymptotic amplitudes of the dark and bright solitons are given respectively by A^i = 1.0 
and A b2 = 0.92 and hence the former is a black soliton. The figure clearly shows the 
solitonic behavior of the solution. The dark soliton suffers a positive phase shift whereas 
the bright soliton suffers a negative phase shift. The formulas Ax\ for the dark soliton 
and Ax 2 for the bright soliton for the phase shifts are given respectively by 

' <n ~ 2 + 2 ' a x = V / K 3 p 2 (l + np 2 ) sin 7l , < 71 < 7r, (4.44a) 



Axi = In 

ai 

Ax 2 = In 

a 2 



sin \ ( 7l - 72 ) 
sin \ (72 + 71) 



sin \ (72 - 71) 



a 2 



a/k 3 p 2 (1 + np 2 ) sin 72, n < 72 < 2tt. (4.446) 



As in the case of the dark-dark solitons, one can see that Axi > and Ax 2 < 0. In the 
present example, Axi = 0.70 and Ax 2 = —0.36. 

4-3. Dark N -soliton solution 

The preceding analysis reveals that the asymptotic form of the iV-soliton solution will be 
represented by a superposition of n dark solitons and N — n bright solitons where n is an 
arbitrary nonnegative integer in the interval < n < N . The derivation of the large time 
asymptotic for the general iV-soliton solution can be done following the similar procedure 
to that used for the two-soliton case. Hence, we outline the result. We address the dark 
soliton solutions satisfying the conditions k > and aj > (j — 1,2, ...,N). The analysis 
for the bright soliton solutions as well as an arbitrary combination of dark and bright 
solitons can be carried out in exactly the same way. 

To begin with, we order the magnitude of the velocity of each soliton as c\ > c 2 > 
... > > 0. We take the limit t — > —00 with 9 n being finite. Since in this limit, \zj\ — > 
for j < n and \zj\ — > 00 for n < j, we find that the leading-order asymptotic of the tau 
function f = J'n from (3.1) with (3.2) can be written in the form 

f N - \(c jk ) n+1 < J>k < N \ fl ( 1 + J^fa^ z n A . (4.45a) 

jtn+l ^ \{Cjk)n+l<j,k<N\ J 

Here, (cjk) is a matrix of Cauchy type given by 

c^ = ^— % l<J,k<N. (4.456) 

pj + pi 
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Referring to the well-known Cauchy's formula, the determinant of the matrix (cjk) is 
evaluated as 

\! \ I TT/ • \ ^lm<j<k<n(Pj ~ Pk)(P*j ~ P*k) 1 ^ , r /. s 

|(cjfc)m<j,fc<n| = ]_]_(« - iPj) =yp 7i | „^ , 1 < m < n < iV. (4.45c) 



If we use (4.45c), we have 



<-m<j,k< 



n(Pj+P*k) 



jk)n<j,k<N | 



K — ipr, 



cxp 



\{Cjk)n+l<j,k<N\ Pn + Pn 

Substitution of (4.46) into (4.45) now gives 



Pn + P*i 



N 



Pn + Pj 



N 

Y In 



. (4.46) 



v 



/iV ~ |(c i fc)„+l<j,fc<Ar| JJ (ZjZ?) ( 1 + 

j=n+l 



K ^Pn / i * 



where 



Z n — Z n ex P 



TV 



E 



j=n+l 



Pn + p„ 

Pn +P*j 
Pn ~ Pj 



(4.47a) 



(4.476) 



The leading-order asymptotic of g n in the limit of t — > — oo can be derived in the same 
way. It takes the form 

9N ~ \(c> k ) n+1 < ]>k < N \ n (l " ^3^ Z 'n Z 'n) > ( 4 - 48 °) 

j=n+1 V Pn^r P n Pn J 



where 



Pj + Pi Pi 



(4.486) 



The asymptotic form of the dark iV-soliton solution follows from (2.1), (4.47) and (4.48). 
It reads 



un ~ pexp {i {kx — out + <p n -*) } 



1 _ Pn I I * 

1 _|_ K~ip„ ii* ' 

1 + Pu+Pn 



with 



0i" 



arg 



N 



n (— ^ 



(N -n)n. 



This expression can be rewritten in terms of the one-soliton solution as 

u N ~ exp Ul (0 n + A0<">), 



(4.49a) 



(4.496) 



(4.50a) 
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with 



N 



A ^r ) = - E ln 

j=n+l 



Pn+P*; 



(4.506) 



Pn ~ Pj 

By a similar asymptotic analysis, we can derive the asymptotic form of un in the limit 
of t — > +00. We find that 



^~exp(i0( + )) Ml (^ + Ae ) ), 



with 



n— 1 

A6>W = -J> 



Pn +pJ 



- Pj 



arg 



'n-l 



+ (n - 1)tt. 



(4.51a) 

(4.516) 
(4.51c) 



n — 

We see from (4.50) and (4.51) that in the rest frame of reference, the asymptotic form of 
the dark iV-soliton solution can be represented by a superposition of N independent dark 
one-soliton solutions, the only difference being the phase shifts of each soliton caused by 
the collisions. It follows from (4.50b) and (4.51b) that the formula for the total phase 
shift of the nth soliton is given by 



Ax - = T ( E ln 



\j=n+l 



Pn + Pj 



Pn - Pj 



n-l 

i=i 



Pn+P*i 



Pn - Pj 



n = 1, 2, N. 



(4.52) 



As in the two-soliton case, we can rewrite the above formula in terms of the variables •jj 
defined by (4.6) and (4.7). Explicitly, 



Ax, r 




sini( 7n + 7j ) 



sin|(7n-7j) 



n-l 



sini(7„ + 7j ; 



sini( 7n - 7j ) 



\/ K 3 p 2 (l + /tp 2 ) sin 



7n 



< 7„ < 7T, 



n 



1,2, ...,7V. 



(4.53) 



The formulas (4.52) and (4.53) reduce to (4.36), (4.37), (4.42) and (4.43) for the special 
case of iV = 2. They clearly show that each soliton has pairwise interactions with other 
solitons, i.e., there are no many-particle collisions among solitons. This feature is common 
to that of the bright iV-soliton solution considered in I. 
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5. Concluding remarks 

In this paper, the system of bilinear equations reduced from the FL equation has been 
derived and used to construct the dark A-soliton solution. The corresponding A-soliton 
solution derived in [7] using the Backhand transformation follows from our solution (2.1) 
with (3.1) and (3.2) if one introduces the angular variables jj according to the relations 
(4.7). We have found that unlike the bright soliton solutions obtained in I, the complex 
amplitude parameters pj are subjected to the constraints (3.2c) which have prevented the 
proof of the solution. To overcome this difficulty, we have employed a trilinear equation 
in place of one of the bilinear equations, in addition to an auxiliary variable r in (3.2c). 
As a byproduct, this trilinear equation has led for the first time to a simple formula for 
the dark A-soliton solution of the derivative NLS equation on the background of a plane 
wave. Note that the dark soliton solutions on a constant background [13-19] stem simply 
from the above-mentioned solution in the zero limit of the wavenumber k. However, this 
limiting procedure is found to be unable to perform for the dark iV-soliton solution of the 
FL equation due to the singularity of the dispersion relation. 

We have seen that the soliton solutions presented here exhibit several new features. 
Specifically, both the dark and bright solitons exist depending on the sign of the wavenum- 
ber k and that of the real part of the complex amplitude parameter. Of particular interest 
is the existence of an algebraic dark soliton which appears only in the case of negative k. 
Finally, the asymptotic analysis of the two- and general A-soliton solutions has clarified 
their structure and dynamics. In particular, the latter solution has been shown to include 
n dark solitons and N — n bright solitons on nonzero background with n being an arbitrary 
nonnegative integer not exceeding N. The application of the results summarized above to 
nonlinear fiber optics will be an interesting issue to be studied in a future research work. 
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